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Abstract. Recent observations give some clues that the 
lenses discovered by the microlensing experiments in the 
direction of the Magellanic Clouds may be located in these 
satellite galaxies. We re-examine the possibility that self- 
lensing alone may account for the optical depth measured 
towards the Large Magellanic Cloud (lmc). We present 
a self-consistent multi-component model of the lmc con- 
sisting of distinct stellar populations, each associated to 
a vertical velocity dispersion ranging from 10 to 60 km/s. 
The present work focuses on showing that such dispersions 
comply with current 20 — 30 km/s limits set by observa- 
tion on specific lmc populations. We also show that this 
model reproduces both the 1 — 2 x lO^*" observed optical 
depth and the event duration distribution. 

Key words: Galaxy: halo, kinematics and dynamics, stel- 
lar content - Cosmology: dark matter, gravitational lens- 
ing 

Several collaborations (Alcock et al. 1993, Aubourg et 
al. 1993) are searching for galactic dark matter through 
the use of gravitational microlensing (Paczyhski 1986) to- 
wards the Magellanic Clouds. Events have been observed, 
for which location and mass cannot be determined inde- 
pendently. The current results do not yet yield a coherent 
explanat ion: half of the hal o of the Milky Way in 0.5 Mq 
objects ( Alcock et al. 1997 ) would require a puzzling star 
formation history, whereas traditional models of the lmc 
do not predict a self-lensing optical depth high enough to 
account for all the observed events (Gould 1995). The only 
events with additional information all seem to be located 
in the Clouds themselves (Bennett et al. 1996, Palanque- 
Delabrouille et al. 1998, Afonso et al. 1999), which makes 
it worthwhile to re-examine the experimental constraints 
on the Clouds kinematics and explore more thoroughly 
models of the lmc. After reviewing the observational con- 
straints on the LMC kinematics (section 1), we show, in 
section 2, the existence of an age bias: the stars used to 
derive these constraints are on average both younger and 



slower than the majority of the lmc objects. Wc then use a 
Monte Carlo simulation to show that a maximum velocity 
dispersion of 60 km/s reproduces the kinematic observa- 
tions (section 3) and the microlensing results (section 4). 

1. Present observational constraints 

The bulk of the mass of the lmc resides in a nearly face- 
on disk, with an inclination usually taken to equal the 



canonical value of « = 33° (Westcrlund 1997), although 
both lower (27°) and higher (up to 45°) values have also 
been derived from morphological or kinematical studies 
of the LMC. This disk is observed to rotate with a circu- 
lar velocity Vc ^ 80 km/s out to at least 8° from the 
LMC center ( |Schommer et al. 1992| ). If all the stars belong 
to the same population, with a vertical {i.e. perpendicu- 
lar to the disk) velocity dispersion aw, the microlensing 
optical depth of such a disk upon its own stars is given 



by T ~ 2cr^ sec i / c ( [Gould 1995)) . Considering the mea- 
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sured velocity of lmc carbon stars (Cowley & Hartwick 
1991), Gould (1995) assumed aw = 20 km/s as a typi- 
cal velocity dispersion for lmc stars. He thus concluded 
that T ^ 10~ ^, i.e. tha t self-lensing (first suggested by 
|Sahu 1994] and |Wu 1994[ ) contributes very little to the ob- 
served optical depth towards this line of sight. 

Carbon stars however may not be the ultimate probe 
to infer the velocity dispersion of LMC populations: they 
actually comprise various ill-defined classes of objects 
(Menessier 1999), and their prevalence is a complex func- 
tion of age, metallicity and probably other factors (Gould 
1999). 

Both observational and theoretical arguments favour 
the existence of a wide range of velocity dispersions among 
the various LMC stellar populations. To commence, Mea- 
theringham et al. (1988) have determined the radial ve- 
locities of a sample of planetary nebulae (FN) in the lmc. 
They measured a velocity dispersion of 19.1 km/s, much 
larger than the value of 5.4 km/s found for the HI. This 
was interpreted as being suggestive of orbital heating and 
diffusion operating in the lmc in the same way as it is ob- 
served in the solar neighbourhood. Then, the observations 
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of Hughes et al. (1991) show clear evidence for an increase 
in the velocity dispersion of long period variables (LPV) 
as a function of their age. For young LPVs, the velocity 
dispersion is 12 km/s whereas for old LPVs, it reaches 35 
km/s. More recently, Zaritsky et al. (1999) found a ve- 
locity dispersion of a = 18.4 ± 1.4 km/s for 190 vertical 
red clump (VRC) star^ whereas for the red clump (RC), 
they measured a value of ct = 32.2 ± 3.8 km/s on a sample 
of 75 objects (throughout this paper, error bars are con- 
verted from Zaritsky's 95 % confidence levels to standard 
la). A general trend appears: the velocity dispersion is an 
increasing function of the age. Just like for our own Milky 
Way, stars of the LMC disk have been continuously un- 
dergoing dynamical scattering by, for instance, molecular 
clouds or other gravitational inhomogeneities. This results 
in an increase of the velocity dispersion of a given stellar 
population with its age, as will be further discussed in sec- 
tion 3. Notice that the main argument in disfavour of a 
LMC self-lensing explanation is precisely the low value of 
the measured vertical velocity dispersions. However, the 
stellar populations so far surveyed predominantly consist 
of red giants. They are shown in the next section not to 
be representative of the bulk of the LMC disk stars, and 
actually biased towards young ages: they are on average 
~ 2 Gyr old, to be compared to an LMC age of ^ 12 Gyr. 

2. The age bias 

The red clump population will illustrate the main thrust 
of our argument. Clump stars have burning helium cores 
whose size is approximately independent of the total mass 
of the object. They also have the same luminosity and 
hence they spend a fixed amount of time ths in the clump, 
irrespective of their mass m. Such objects are evolved 
post-MS stars, which does not mean that they are neces- 
sarily old. We have assumed a Salpeter Initial Mass Func- 
tion for the various LMC stellar populations 



dm 



-(1 + a) 



(1) 



with a = 1.35. The stellar formation history has been 
borrowed from Geha et al. (1998). Their preferred model 
(e) corresponds to a stellar formation rate J-{t) that has 
remained constant for 10 Gyr since the formation of the 
LMC 12 Gyr ago. Then, two Gyr ago, T{t) has increased by 
a factor of three. The number of stars that formed at time 
t and whose mass is comprised between m and m + dm 
may be expressed as 



d^N 
dm dt 



(2) 



We have assumed a mass-luminosity relation L oc m^ on 
the MS so that the stellar lifetime may be expressed as 



^ see Zaritsky et al. (1999) and Beaulieu and Sackett (1998) 
for a definition of RC and VRC stars. 



tms("i) = 12 Gyr/? 



(since r oc m/L). With these 



oversimplified but natural assumptions, a star whose ini- 
tial mass is < 1 Mq is still today on the MS and cannot 
have reached the clump. Conversely, a heavier star with 
m > 1 Mq may well be today in a helium core burning 
stage provided that its formation epoch lies in the range 
between t — —TMsim) (the object has just begun core 
helium burning) and t = — tms(w) — rHe(w) (the star is 
about to leave the red clump). The number of RC stars 
observed today with progenitor mass in the range between 
m and m -I- dm is therefore given by 



dVRc = ^(-TMs(™)) X TO (1 + ") dm X thc • 



(3) 



To get more insight into the age bias at stake, we can 
parameterize the progenitor mass m in terms of the age 
T = tms(to). The previous relation simplifies into 



RC 



dr 



•^(-t) THc 



^(7-1) 



(4) 



where "f = a/ (/3— 1). This may be directly compared to the 
age distribution of the bulk of the LMC stars that goes like 
!F{—t). With a Salpeter mass function and (3 = 4.5, we 
get a value of 7 = 0.4. The excess of young RC stars goes 
as and the bias is obvious. Other IMF are possible 

and a spectral index as large asa^^/? — 1~3.5 would 
be required to invalidate the effect. HST data analyzed 
by Holtzman et al. (1997) nevertheless point towards a 
spectral index a that extends from 0.6 up to 2.1 for stars 
in the mass range 0.6 < m < 3 Mq. The average value 
corresponds actually to a Salpeter law. 

There has been furthermore a recent burst in the LMC 
stellar formation rate. In order to model it, we may express 
the total number of today's RC stars as an integral where 
the progenitor mass m runs from mi = 1 M© up to the 
tip of the IMF whose actual value is irrelevant and has 
been set equal to infinity here for simplicity. Notice that 
the specific progenitor mass TO2 ~ 1.7 Mq corresponds 
to stars born 2 Gyr ago, when the stellar formation rate 
increased by a factor of 3. Stars which formed before that 
epoch will be referred to as old. Their number is given by 



ATOld 



(5) 



On the other hand, the number N^q^^ of young clump 
stars is obtained similarly, with masses in excess of TO2. 
We readily infer a fraction of young stars 



V, 



younji 
RC 



/V] 



RC 



2 -I- (TO2/mi)" 



0.751 . 



(6) 



Three quarters of the clump stars observed today in the 
LMC have thus formed less than 2 Gyr ago, during the 
recent period of stellar formation mentioned above. Inte- 
grating tms over the RC population 



(r) 



1 



Tms dN^c 



(7) 
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yields the average age 
(r) = (12 Gyr) x 



a 



2 m. 



a + 1 



2 mn 



(8) 



This gives a numerical value of ~ 1.95 Gyr. We thus 
conclude that today's clump stars are, on average, much 
younger than the lmc disk. 

3. Distributions of velocity dispersions 

This simple analytical result has been checked by means of 
a Monte Carlo study. We have randomly generated a sam- 
ple of 10^ LMC stars. The progenitor mass was drawn in 
the range 0.1 < m < 10 Mq according to a Salpeter law. 
The age of formation was drawn in the range —12 Gyr 
< i < according to the stellar formation history J-{t) 
favoured by Geha et al. (1998). The vertical velocity dis- 
persion aw was then evolved in time from formation up 
to now according to Wielen's (1977) relation: 



(9) 



This purely diffusive relation is known to be inadequate to 
describe velocity dispersions in our Galaxy (Edvardsson et 
al. 1993). We will however use it in our model, as heat- 
ing processes in the lmc may be different than those in 
the galaxy. The lmc is indeed su bject to tidal heating by 
the Milky Way ( Weinberg 1999| ) and has most probably 
suffered encounters with the SMC. Although this simple 
relation lacks a theoretical motivation, it will be shown to 
account for several features of the velocity distributions in 
the LMC, without being at variance with any observation. 
The initial velocity dispersion ao was taken to be 10 km/s, 
and the diffusion coefficient in velocity space along the ver- 
tical direction Cw to be 300 km^ s^^ Gy^^ so that our 
oldest stars have a vertical velocity dispersion reaching up 
to o■^A'^ — 60 km/s. For each star, the actual vertical 
velocity was then randomly drawn, assuming a Gaussian 
distribution with width aw- 

In order to compare our Monte Carlo results with the 
Zaritsky et al. (1999) measurements of the radial veloci- 
ties of LMC clump stars, we selected two groups of stars 
according to their position in the HR diagram. Following 
Zaritsky et al., we use their colour index 



C = 0.565 (5-/) + 0.825(C/- F + 1.15) 



(10) 



so that the RC population is defined by 3.1 < C < 3.4 
with a magnitude 19 < < 19.3 whereas the VRC 
stars have the same colour index C and brighter mag- 
nitudes 18 < V < 18.75. In order to infer the colours and 
magnitudes of the stars that we generated, we used the 
isochrones computed by Bertelli et al. (1994) for a typical 
LMC metallicity and helium abundance oi Z = 0.008 and 
Y = 0.25. 

A random sample of 190 stars that passed the VRC 
selection criteria is presented in Fig. Q where the vertical 



Vertical Red Clump 




RMS = 1 8 km/s 



Velocity (km/s) 

Fig. 1. Velocity distribution for a sample of 190 vertical 
red clump stars that have been generated by the Monte 
Carlo discussed in the text. That histogram is similar to 
Fig. 10 of Zaritsky et al. (1999). A velocity dispersion of 
18 km/s is found for the full sample (solid smooth curve). 



velocities are displayed. This histogram may be compared 
to Fig. 10 of Zaritsky et al. (1999) where no VRC star is 
found with a velocity in excess of 60 km/s. With the full 
statistics, our Monte Carlo generated a population of ~ 
2,900 VRC objects whose vertical velocity distribution has 
a RMS of ~ 18 km/s. The agreement between the Zaritsky 
et al. observations and our Monte Carlo results is notewor- 
thy. The average age of our VRC sample is ~ 0.87 Gyr. 



Red Clump 




RMS . 23 km/s 



100 150 

Velocity (Itm/s) 



Fig. 2. Like in the previous figure, a distribution of 75 
red clump stars is now featured. We inferred a velocity 
dispersion of 23 km/s for the full sample (solid smooth 
curve). Our distribution is similar to that presented in 
Fig. 11 of Zaritsky et al. (1999). No star exhibits a velocity 
larger than 70 km/s. 



We also selected a random sample of 75 RC stars whose 
velocity distribution is featured in Fig. |^. Even with a dif- 
fusion coefficient as large as Cw = 300 km^ s^^ Gy^^ so 
as to comply with a large lmc self-lensing optical depth, 
our full statistics of 18,000 RC objects has a velocity dis- 
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persion of ~ 23 km/s. This is sliglitly below the value of 
a = 32.2±3.8 km/s quoted by Zaritsky et al. Observations 
are nevertheless fairly scarce with only 75 RC stars. When 
Zaristsky et al. fitted a Gaussian to the RC radial veloc- 
ity distribution featured in the Fig. 11 of their paper, they 
obtained a 95 % C.L. dispersion of cr = 32^i[g km/s with 
a large uncertainty. Our Monte Carlo velocity dispersion 
of 23 km/s is definitely compatible with that result. We 
infer an average age for the RC population of ^ 1.8 Gyr 
to be compared to our analytical result of ^ 1.95 Gyr. 
This agrees well with Beaulieu and Sackett's conclusion 
that isochrones younger than 2.5 Gyr are necessary to fit 
the red clump. Notice finally that our age estimates for 
these various clump populations are in no way related to 
LMC kinematics. They merely result from the postulated 
Salpeter IMF, the Geha et al. preferred stellar formation 
history and the Bertelli et al. isochrones. 

With this model, 70% in mass of the lmc disk consists 
of objects whose vertical velocity dispersion is in excess of 
25 km/s, although the average vertical velocity dispersion 
of RC stars, for instance, is only ~ 23 km/s. 

What about the other measurements? The velocity 
dispersion of PNs has been found equal to 19.1 km/s 
(Meatheringham et al. 1988). These authors estimate that 
the bulk of the PNs have an age near 3.5 Gyr. They also 
note that younger objects are present down to an age of 
order 0.5 — 1.3 Gyr. Meatheringham et al. come finally 
to the conclusion that the indicative age of the FN popu- 
lation is 2.1 Gyr. This value agrees well once again with 
our analytical estimate. Our Monte Carlo gives a slightly 
larger value of 2.4 Gyr for the age of the PNs, with a veloc- 
ity dispersion of 24.7 km/s. Because the observed sample 
contains 94 objects, the measured value of 19.1 km/s suf- 
fers presumably from significant uncertainties. 

Quite interesting also are the measurements by Hughes 
et al. (1991) of the velocity dispersions of LPVs as a func- 
tion of their age. Their sample of 63 "old" LPVs has a 
velocity dispersion of cr = 35l4° km/s. For the bulk of 
the LMC populations, we obtain an average velocity dis- 
persion of ~ 37 km/s. The problem at stake is actually 
the age of those old LPVs. These stars indeed display an 
age-period relation. However, Hughes et al. derived this re- 
lation from kinematics considerations, using precisely Eq. 
^, and postulating the same diffusion coefficient as in the 
Milky Way. They thus inferred an average age of 9.5 Gyr. 
Finding instead the position of these stars in a colour- 
magnitude diagram and using LMC isochrones would have 
led to a clean determination of the age-period relation. A 
direct determination of the age of LPVs is nevertheless 
spoilt by a few biases. Some LPVs are carbon stars and 
the ejected material around them may considerably dim 
their luminosities. These stars may also pulsate on an har- 
monic of the fundamental mode. Both effects lead to an 
under-determination of their luminosity and hence to an 
overestimate of their age (Menessier 1999). As a matter 
of fact, Groenewegen and de Jong (1994) conclude that 



LMC stars whose progenitor mass is less than 1.15 M© 
never reach the instability strip on the AGB. This yields 
an upper limit on the age of LPVs of ^ 7.3 Gyr, in clear 
contradiction with the average age of 9.5 Gyr inferred by 
Hughes et al. for old LPVs. 

Finally, Schommer et al. (1992) have obtained a ve- 
locity dispersion of 21 — 24 km/s for 9 old LMC clusters. 
Their large la error of ^ 10 km/s is due to the small 
size of the sample. It is not clear whether or not these 
clusters have formed in the disk. If they nevertheless had, 
they would have undergone a fairly restricted orbital heat- 
ing with respect to the LMC stars. Those systems and the 
giant molecular clouds have actually comparable masses 
and the energy exchange between them does not result 
in a significant increase of the velocity dispersion of the 
clusters unlike what happens to the stars. 

4. Multi-component model of the LMC 

We model the LMC to contain several stellar populations, 
each associated with a different velocity dispersion aw,i 
which has evolved according to Eq. |[ 

We describe each of the ten components of our model 
by an ellipsoidal density profile 



Pi{R,z) 



A,, 



i?2 + z'/{l~ei) 



(11) 



up to a cut-off radius -Rmax = 15 kpc (Aubourg et al. 
1999). The multi-component model based on these pro- 
files is self-consistent in the sense that it satisfies Pois- 
son equation and results in a flat rotation curve with the 
desired Vc of 80 km/s. We define the set of aws so as 
to sample linearly the range between ctq = 10 km/s and 
60 km/s (see previous section). The parameters 



^MAX 



Ki and the ellipticities are determined so that the model 
reproduces the set of velocity dispersions uw,i and surface 
mass densities where dTn/dai oc aiJ-{t) with J^{t) the 
stellar formation history of the LMC mentioned in section 
2. Assuming a typical M/L of 3, which is a free parameter 
in our model, we reproduce the observed surface bright- 
ness of the LMC. 

For a given distribution of objects, one can com- 
pute the total self-lensing optical depth r and the event 
rate F. Both quantities are integrated on all deflectors 
and sources, considering that only main sequence stars 
brighter than V = 2Q and red giants can be potential 
sources, since they are the only objects bright enough to 
be visible in microlensing surveys. The computation of F 
requires an estimate of the relative transverse velocity of 
deflector and source, for which we have assumed an hor- 
izontal velocity dispersion equal to the vertical one pre- 
dicted by the model. Details of this computation can be 



found in (Aubourg et al. 1999). 

For the model described above, one obtains r = 
9.3 X 10"^ and F = 3.5 x 10"'^ yr"!. This can be com- 
pared to the EROS and MACHO optical depths, respectively 
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8.2 X 10-8 dAnsari ct al. 1996D and 2.9tli ^ 10"'^ (Alcock 
et al. 1997). A combination of those two results yields an 



+ 1.3 

-0.8 



X 



10"'^ ( [Bennett 199^ , 



average optical depth of 2.1 
but preliminary MACHO results from their five-year anal- 
ysis (Sutherland 1999) hint to a reduced optical depth as 
compared to their two-year analysis. The model predic- 
tion is thus in good agreement with the results obtained 
so far from microlensing experiments. 

Another relevant prediction of the model is the dis- 
tribution of event durations, dT/dAt. Figure |3| illustrates 
this prediction for our model, along with the distribution 
of observed macho events. 




Cumulative duration distribution 

Distribution from MACHO events 

Distribution from model 
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Fig. 3. Predicted distribution of event durations dT/dAt, 
superimposed with the MACHO experimental distribution. 
The events are those presented by MACHO at the IVth Mi- 



crolensing Workshop (Cook 1998), corrected for blending 
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and efficiency using the formulae in Alcock et al. 1997 



Our model thus reproduces both the total observed 
optical depth towards the LMC and the observed event 
duration distribution, while complying with the velocity 
dispersion measurements. A self-lensing interpretation of 
all the microlensing events observed so far towards the 
LMC thus appears to be a plausible explanation. 
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